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In order to clarify possible multipole states of Sm-based filled skutterudite compounds, we investigate 
multipole susceptibility of a multiorbital Anderson model dynamically coupled with Jahn-Teller phonons 
by using a numerical renormalization group method. Here we take a procedure to maximize the multipole 
susceptibility matrix to determine the multipole state. When the electron-phonon coupling term is simply 
ignored, it is found that the dominant multipole state is characterized by 2u octupole for the F^^ quartet 
ground state, while the low-temperature phase is governed by magnetic fluctuations for the F J doublet 
ground state. When we include the coupling between / electrons in degenerate F^^ (Cu) orbitals and Jahn- 
Teller phonons with Eg symmetry, the mixed multipole state with 4u magnetic and 5u octupole moments 
is found to be dominant with significant 3g quadrupole fluctuations at low temperatures. 
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I. Introduction 

Recently, /-electron compounds including plural number 
of / electrons per ion have attracted renewed attention due 
to emergence of exotic magnetism and unconventional su- 
perconductivity in the research field of condensed matter 
physics.' In particular, filled skutterudite compounds, ex- 
pressed as RT4X12 with rare-earth atom R, transition metal 
atom T, and pnictogen X, have been focused, since this ma- 
terial group provides us a platform for systematic research 
of magnetism and superconductivity of /"-electron systems 
with n>2,^'^ where n denotes the number of / electrons. 

Among various kinds of RT4X12, interesting results have 
been reported for Sm-based filled skutterudites. SmRu4Pi2 
is known to exhibit a metal-insulator (MI) transition at 
7mi=16.5 K,** which is confirmed to be second order from 
the specific heat measurement.^ Further works have revealed 
that the MI transition occurs in two steps. Namely, we find 
three states in the phase diagram for SmRu4Pi2 in the plane 
of temperature and magnetic field. ^ The phase I (T>Tmi) is 
paramagnetic, while the phase III is considered to be antifer- 
romagnetic ordered state from NMR measurements.^-'^ 

Concerning the phase II between the phases I and III, an 
intriguing possibility of octupole ordering has been proposed 
by Yoshizawa et al. from the elastic constant measurement, 
showing the breaking of time reversal symmetry in the phase 

II. '" Recent muon spin relaxation measurement is consistent 
with this octupole scenario, since the static internal field has 
been found to grow at Tmj." NMR measurement of ■^'^P has 
also suggested that the static internal field occurs below Tmi 
and the anomalous variation of the internal field is considered 
to be the evidence of the ordering of octupole.^- '^ 

In order to understand why such exotic octupole ordering 
appears, it is important to clarify the CEF ground state. From 
the specific heat measurement, for SmRu4Pi2 and SmOs4Pi2, 
the CEF ground state is F^y quartet,'-' consistent with theoret- 
ical calculation of the CEF energy level of Th group. Note 
that in Th group is equal to F^ in Oh - However, the CEF 
energy scheme seems to be changed even among the same 
Sm-based filled skutterudites. In fact, for SmFe4Pi2, F^^ dou- 
blet ground state has been suggested. '^''^ For SmOs4Sbi2, it 
has been found to be Fgy quartet in Ref. 17, while the CEF 



ground state was considered to be F^ doublet in Ref. 18. 
Quite recently, magnetization measurement has been per- 
formed and the observed anisotropy seems to confirm the F^^. 
quartet ground state in SmOs4Sbi2.'^ 

Another interesting topic on Sm-based filled skutteru- 
dites is magnetically robust heavy-fermion behavior in 
SmOs4Sbi2.'^ It has been suggested that this peculiar phe- 
nomenon is due to the non-magnetic Kondo effect originat- 
ing from phononic degrees of freedom.^" Since rare-earth ion 
is suiTounded by the cage composed of twelve pnictogens in 
filled skutterudites, there is a possibility that rare-earth ion 
moves around potential minima in off-center positions inside 
the pnictogen cage. This is called the rattling, which is con- 
sidered to be one of important ingredients with significant in- 
fluence on electronic properties of filled skutterudites. As a 
natural extension of two-level Kondo problem,^'- Hattori et 
al. have analyzed four- and six-level Kondo models. '^^ Their 
results seem to explain the magnetically robust heavy-fermion 
behavior in SmOs4Sbi2. Concerning the positions of potential 
minima, neutron scattering experiment for PrOs4Sbi2 has in- 
dicated that the charge distribution of Pr ion extends in the 
[1,1,1] direction at high temperatures,^^ suggesting that the 
eight-level Kondo model should be discussed. 

Regarding the problem of rattling, there is another impor- 
tant aspect of geometrical degree of freedom. From the mea- 
surement of elastic constant, it has been suggested that the off- 
center motion has degenerate symmetry. Inspired by this 
suggestion, the present author has considered a multiorbital 
Anderson model with a linear coupling between degenerate /- 
electron orbitals with eu symmetry and dynamical Jahn-TeUer 
phonons with Eg symmetry. Numerical analysis of this 
model for n=2 has revealed that quasi-Kondo behavior occurs 
due to the release of an entropy log 2 of the vibronic ground 
state, originating from the clockwise and anti-clockwise rota- 
tional modes of dynamical Jahn-Teller phonons.^** Such geo- 
metrical degree of freedom does not depend on the detail of 
potential minima, since only the rotational direction becomes 
important. It is interesting to pursue a possibility of Kondo- 
like behavior due to the dynamical Jahn-Teller effect in filled 
skutterudites. 

In this paper, we discuss multipole properties of Sm-based 
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filled skutterudites on the basis of the multiorbital Anderson 
model. First we introduce the effective local /-electron model 
obtained by the expansion in terms of 1/A based on a j-j cou- 
pling scheme, where A is the spin-orbit coupling. It is found 
that the ground state is changed between doublet and 
quartet, depending on Coulomb interaction and/or spin- 
orbit coupling. This fact is consistent with a possible change 
of the CEF ground state in SmT4Xi2. Then, we numerically 
analyze the multiorbital Anderson model. When the ground 
state is F^^ quartet, the dominant multipole state is charac- 
terized by 2u octupole. For the F^ doublet ground state, the 
low-temperature phase is governed by magnetic fluctuations. 
When we consider the coupling with Jahn-Teller phonons, the 
mixed multipole state with 4u magnetic and 5u octupole mo- 
ments becomes dominant at low temperatures and quadrupole 
fluctuations are also significant. 

The organization of this paper is as follows. In §2, we intro- 
duce the effective model which describes well the low-energy 
local /-electron states. On the basis of the model, we dis- 
cuss the CEF energy levels of Sm^+ from the viewpoint of 
the competition between Coulomb interaction and spin-orbit 
coupling. In §3, the multiorbital Anderson model is defined 
and we explain the quantities to be measured and the method 
of calculation. In §4, we show numerical results for the cases 
without and with Jahn-Teller phonons. In particular, the dom- 
inant multipole states are explained in detail. Finally, in §5, 
we discuss relevance of the present results to actual Sm-based 
filled skutterudites. Throughout this paper, we use such units 
as h=kQ=l and the energy unit is set as eV. 

2. Local /-electron State 

2. 7 Local effective model 

The local /-electron Hamiltonian is given by'"' 

-ffloc = Hso + Hint + HqeF, (1) 

where Hgo is the spin-orbit coupling term, given by 

-^so ~ A ^ ^ Cmja,m' ,a* frnaf^T^'o'' ' (2) 

Here A is the spin-orbit coupling and f^u denotes the anni- 
hilation operator for / electron with spin a and angular mo- 
mentum m(=— 3,- • • ,3) and a=+l (—1) for up (down) spin. 
The matrix elements are expressed by 

Cm,(T,m,cr ~ VflG 

Cm+o-,-a,m,cr = -\/l2 - TO(m + cr)/2, 

and zero for other cases. 

The second term H-^nx indicates the Coulomb interactions 
among / electrons, expressed by 

-^m 1 , m2 ,7713 ) 7714 /m 1 (71 /rn2 ctq. f'ms cr2 /77i4 (Ti •> (4) 

where the Coulomb integral Imi,m2,7n3,m4 is given by 

6 

-^mi ,m2,7n3 ,m4 — 

fe=o 

Note that the sum is limited by the Wigner-Eckart theorem 
to even values (fc=0, 2, 4, and 6). Here F'' is the radial inte- 
gral for the fc-th partial wave, called Slater integral or Slater- 
Condon parameter^**' and Cfe is the Gaunt coefficient.^'^' 



The Slater-Condon parameters and the spin-orbit interaction 
are determined so as to reproduce the spectra of rare-earth 
and actinide ions for each value of n, but the Slater-Condon 
parameters are considered to be distributed among 1 eV and 
10 eV, while A is in the order of 0.1 eV. In this paper, for 
convenience, we parameterize F'^ as 

F° = IQU, = 5U, F^ = 3U, F^ = U, (6) 

where U denotes the order of the Hund's rule interaction 
among / orbitals, which is in the order of a few eV. 
The CEF term Hqef is given by 

HcEF = ^ Bm,m'fmafm'a, (7) 

where Bm,m' is determined from the table of Hutchings for 
angular momentum J=£=3,^^ since we are now considering 
the potential for / electron. For filled skutterudites with Th 
symmetry,'^ Bm,m' is expressed by using three CEF parame- 
ters Bl S§, and ^l as 

B3,3 = B-3,-3 = 1805° + 1805°, 
B2,2 = B_2.-2 = -420B2 - lOSOBg, 
= = 605^ + 27005°, 

Bo,o = 360B2 - aeooBg, 

Ba'-i = B_3,i = Q0VT5{B2 - 21BI), (8) 
B2,-2 = 300S2 + 75605°, 
B3,i= B-3,-i = 24VTEbI, 
B2.0 = S_2.o = -48VT5B|, 
Bi,-i = -B3,-3 = 360B|. 

Note the relation of Bm,m'=Bm',m- Following the traditional 
notation, we define 

B2 = Wx/F{4), 

B°^ = W{l-\x\)/F{6), (9) 
Bl = Wy/F\6), 

where x and y specify the CEF scheme for Th point group, 
while W determines an energy scale for the CEF potential. 
Concerning F(4), F(6), and F*(6), we choose F(4)=15, 
F(6)=180, and F*(6)=24 for ^=3.34.35 in actual /-electron 
materials, the magnitude of the CEF potential is considered 
to be 10~^~10~'^eV. For the calculation of the CEF energy 
level in this paper, we set W^=— 6x 10~^eV. 

Although H\oc provides us correct CEF ground states irre- 
spective of the number of / electrons and other parameters, 
it is difficult to consider multipole properties and the effect 
of Jahn-Teller phonons. In order to make further steps, it is 
convenient to employ the j-j coupling scheme.^^ Namely, we 
accommodate / electrons in the j=5/2 sextet by considering 
the CEF potential and Coulomb interactions. An advantage 
of the j-j couphng scheme is that we can take into account 
many-body effects using standard quantum-field theoretical 
techniques, since individual /-electron states are clearly de- 
fined. Then, we have developed microscopic theories for mag- 
netism and superconductivity on the basis of the j-j coupling 
scheme by focusing on the potential roles of /-electron or- 
bitals. ^^^^ In particular, it is also possible to develop micro- 
scopic theory for multipole ordering and fluctuations.*^^^ 

However, since the sixth order CEF terms cannot be in- 
cluded in the j=5/2 sextet due to the symmetry reason, con- 
tributions of -Bg and Bq are dropped in the simple /-/ cou- 
phng scheme. It has been a disadvantage of the /-/ coupling 
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scheme that the effect of which is characteristic of Th 
symmetry, cannot be included. In order to improve the situa- 
tion, recently, the modified j-j coupling has been developed 
so as to include the effect of -Bg and i?| terms in the order of 
1/ A due to the degenerate perturbation theory. •^^ The modified 
j-j coupling scheme can correctly reproduce the CEF energy 
levels for n>2 in the realistic intermediate coupUng region 
with X/U in the order of 0.1. 

In order to obtain the effective model H^g, we treat 
^^CEF + ^int as a perturbation to Hso- The details are ex- 
plained in Ref. 29. Here we show i?eff as 

^eS = B^^iifj^f,^ + ^A»u"2,/:i3jii4/^i/^2//^3/jii4) (10) 

where is the annihilation operator for / electron with an- 
gular momentum /x(=— 5/2,- • • ,5/2) in the j=5/2 sextet. The 
modified CEF potential is expressed as 



B = + b(i) 



(11) 



where Bjj'}, denotes the CEF potential for J=5/2 and B^^}, is 
the correction in the order of W^/X. The effective interaction 
in eq. (10) is given by 

7 =7(0) .7(1) (12) 

where -?iil/i2,/*3,/i4 is expressed by three Racah parameters, 
Eq, El, and £'2, which are related to the Slater-Condon pa- 
rameters as 

Eo = F°- (80/1225)i^2 _ (12/441)F'^, 

El = (120/1225)^2 + (18/441)^4, (13) 

E2 = (12/1225)^2 - (1/441)^4. 

Explicit expressions of /'^"^ are shown in Ref. 36. On the other 
hand, iIh,h2,h3,h4, is the correction term in the order of 1/A. 

Here, four comments are in order, (i) Effects of Bg and B| 
are included as two-body potentials in (ii) The lowest- 
order energy of /(I) is \W\U/X. (iii) The effective model 
is vaUd for /-electron compounds with the hybridization 
V smaller than A. (iv) The parameter space in which Hes 
works is determined by the conditions for the weak CEF, 
i.e., \W\/U<^1 and \W\U/X<i:E2. Concerning the value of 
A, we obtain |W|/A<0.02 from eqs. (6) and (13). Note that 
E2 plays a role of the effective Hund's rule interaction in the 
j-j coupling scheme, estimated as JS2~J'h/49.2' Thus, it is 
allowed to use Hcff even for A in the order of 0. 1 eV, when 
\W\ is set as a realistic value in the order of 10""^ eV for ac- 
tual /-electron materials. 

2. 2 CEF energy schemes 

In order to determine the CEF parameters of filled skut- 
terudite compounds, first let us consider the CEF energy level 
for the case of n=2. A typical material is PrOs4Sbi2. From 
experimental results on specific heat,^^ magnetization,^^ and 
neutron scattering, it has been confirmed that the ground 
state of PrOs4Sbi2 is singlet, while the first excited state 
is r4 ^ triplet with very small excitation energy (~10K). 

In Fig. 1(a), we show the CEF energy levels of i?eff for 
the case of n=2. Here we set A=0.1, U=l, W=—6xl0~^, and 
y=0.3, which are considered to be realistic values for Pr-based 
filled skutterudites. Here the minus sign is added in so 
as to be consistent with the appearance of singlet ground 
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Fig. 1. CEF energy levels of //off for (a) n=2 and l/=l, (b) n=3 and C/=l, 
(c) n=5 and (7=1, and (d) n=5 and i7=2.5. We set A=0.1,W=-6xlO-'', 
and y=0.3 for (a)-(d). (e) Phase diagram of the CEF ground state of i/efi 
for n=5, W=—() x 10"'', x=QA, and 2/=0.3. (f) Electron configurations for 
n=\ and 5 in the j-j coupling scheme for B4<0. 



State. It should be noted that two triplet states in Oh, and 
r^, are mixed as r4 and r4 in Th. As understood from 
Fig. 1(a), a;=0.4 seems to correspond to the actual situation of 
PrOs4Sbi2, since the CEF excitation energy is 1 meV, which 
is consistent with experimental results. 

Now we substitute the rare-earth ion from Pr'^+ (n=2) to 
Nd^+ (n=T>) and Sm'^+ (n=5). Here we note that in princi- 
ple, it is not necessary to change the CEF parameters even 
when rare-earth ion is substituted, since the CEF potential is 
given by the sum of electrostatic potentials from hgand ions. 
Then, it is reasonable to use the same CEF parameters also 
for NdOs4Sbi2 and SmOs4Sbi2. In actuality, the electrostatic 
potential may be changed depending on the ion radius of R 
and/or the kinds of T and X. In this sense, x, y, and the abso- 
lute value of |T4^| may be changed among filled skutterudite 
compounds, but there is no reason to change the sign of W. 
Throughout this paper, we keep W as negative. 

Let us discuss the CEF states for n=3.^^ In Fig. 1(b), we 
show the CEF energy levels of H\oc for n=3 by using the same 
parameters as in Fig. 1(a). At x=0.4, the ground state is Fg^ 
quartet and the excited state is another F^-. quartet, consis- 
tent with the experimental result for NdOs4Sbi2.^^ As for the 
CEF excitation energy, it is 0.02eV in the present calculation, 
while it has been experimentally estimated as 220K.^** We can 
obtain the CEF energy for NdOs4Sbi2 consistent with the ex- 
perimental result by using the CEF parameters of PrOs4Sbi2. 

In Fig. 1(c), we show the CEF states for n=5 by using the 
same parameters as in Fig. 1(a). Since we obtain two mul- 
tiplets, and Fg^, for n=5, it is convenient to draw the 
excitation energy. At a;=0.4, the ground state is Fg doublet. 
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while the experimental result for SmOs4Sbi2 has suggested 
Fgy quartet ground state. Here we note that the values of 
U and A can be different for different rare-earth ions. In the 
cases of n=2 and n=3, even if U and A are changed within 
the realistic range, the CEF energy scheme does not change 
qualitatively. However, this is not the case for n=5. 

In Fig. 1(d), we show the results of the CEF energy levels 
for [7=2.5 and A=0. 1 by using the same CEF parameters as 
in Fig. 1(c). It is observed that the x dependence is reversed 
from that in Fig. 1(c), suggesting that the CEF ground state is 
changed, even though we do not change the CEF parameters. 
At x=OA, the ground state is Fg^ quartet and the excitation 
energy is about 3 meV in the present calculation, which is 
comparable with the experimental value of 19K." 

We note that the conversion of the ground state by the 
change of U (or A) is not due to the approximation of i?eff, 
but an intrinsic feature of / electrons, since we could find 
the same results when we diagonaUze simultaneously CEF 
potential. Coulomb interaction, and spin-orbit coupling. By 
changing U and A, we depict the phase diagram for the CEF 
ground state for the case of n=5 and x=OA, as shown in 
Fig. 1(e). Note that in the region of very small A, HeS does 
not work, but it is still available for A:»0.015 for the present 
parameter choice. Roughly speaking, for large A and small U, 
the ground state is Fg^ quartet, while for small A and large U, 
Fg doublet is the ground state. 

In order to understand the change of the CEF ground state, 
it is useful to consider the two limiting situations, A<CC/ and 
X^U, corresponding to the LS and j-j coupling schemes, 
respectively. First we consider the case of the LS coupling 
scheme. The relevant CEF parameter is B^, which is ex- 
pressed as B^=Ail3j{n){r'^), where is a parameter de- 
pending on materials, /3j(n) is the so-called Stevens factor, 
appearing in the method of operator equivalent,^^ and (r^) 
indicates the radial average concerning the local /-electron 
wavefunction. When we assume that A^ is invariant in the 
same material group and (r"*) does not depend on the /- 
electron number, the CEF parameter is basically determined 
by /3j(n), which depend on n, U, and A.^' 

In the LS coupling scheme, /35/2(l)=(l l/7)/3£ and 
/35/2(5)=(13/21)/3^, where /3^=/33(l)=2/495 for £=3. Note that 
/^s /2 {n) is different between the cases of n=l and 5 even for 
the same J=5/2, since L=3 and 5=1/2 for n=l, while _L=5 and 
5=5/2 for n=5, where L and S are total angular momentum 
and spin momentum, respectively. Namely, in the LS cou- 
pling scheme, should have the same sign for r?,= l and 5. 
The CEF ground state for n=l is Fg doublet for iJ4>0, while 
it is quartet for B^<Q. Thus, in the LS coupling scheme, 
it is easy to understand that the ground state for n=5 is the Fg-. 
quartet in the region of a;>0 for the present parameterization 
of Bl=Wx/l^ with negative W. 

Next let us consider another limiting case corresponding 
to the j-j coupling scheme, in which we simply accommo- 
date electrons in the one-electron potentials,^^ as shown in 
Fig. 1(f). For x> with W<Q, i.e., negative B^, we put 
five electrons in the level scheme of lower Fg^ and higher 
F^. Thus, we immediately understand that the ground state 
becomes F^ doublet for n=5 in the j-j coupling scheme, 
when Fg7 quartet is the ground state for n=l. In fact, in 
the j-j coupling scheme, we obtain /35/2(l)=(ll/7)/?^ and 
/35/2(5)=(-ll/7)/3^. The sign of 3% is different between the 



cases of n=l and 5 in the j-j coupling scheme. When we 
change the values of U and/or A from the LS to j-j cou- 
pling regions, we can find the smooth change of (3^ /2 (5) from 
(13/21)/?^ to i-nn)l3e, as found in Ref. 29 

The present result suggests a "non-CEF" scenario to un- 
derstand the change of the CEF ground state in Sm-based 
filled skutterudites. As mentioned in the introduction, for 
SmRu4Pi2 and SmOs4Pi2, the CEF ground state is Fg^ quar- 
tet, while the ground state of SmFe4Pi2 is doublet. In or- 
der to explain the change of the ground state within the LS 
coupUng scheme, it is necessary to change the sign of B^. 
The magnitude of the CEF potential may be changed due to 
the difference in the ion radius of R and the substitution of 
T and/or X in RT4X12, but it seems difficult to imagine that 
the sign of the CEF potential is changed among SmT4Xi2 . As 
understood from Fig. 1(e), it is observed that the ground state 
conversion occurs around the region for realistic values of A 
and U. Namely, due to the change in U and/or A, the ground 
state is easily interchanged for actual Sm materials. 

3. Model and Method 

3. 1 Multiorbital Anderson Hamiltonian 

From the band- structure calculations, it has been revealed 
that the main conduction band of filled skutterudites is Ou with 
xyz symmetry,^" which is hybridized with / electrons in the 
F^ state with symmetry. In order to specify the /-electron 
state, we introduce "orbital" index which distinguishes three 
kinds of the Kramers doublets, two Fg^ and one F^ . Here "a" 
and "b" denote the two Fgy's and "c" indicates the F^. 

Then, the multiorbital Anderson model is given by 

H = Y.^kci,Cka + Y.^Vclj,„ + h.C.) 

ka ka (14) 

where £fc is the dispersion of conduction electrons with F^ 
symmetry, /^^ is the annihilation operator of / electrons on 
the impurity site with pseudospin a and orbital 7, Cfc^ is the 
annihilation operator for conduction electrons with momen- 
tum fc and pseudo-spin a, and V is the hybridization between 
conduction and / electrons with a-^ symmetry. Throughout 
this paper, we set V=Q.Q5 eV. Note that the energy unit of H 
is half of the bandwidth of the conduction band, which is con- 
sidered to be in the order of 1 eV, since the bandwidth has 
been typically estimated as 2.7 eV for PrRu4Pi2.*' Thus, the 
energy unit of H is taken as eV. To set the local /-electron 
number as n=5, we adjust the /-electron chemical potential. 

The last term in eq. (14) denotes the electron-phonon cou- 
pling. Here, the effect of rattling is included as relative 
vibration of surrounding atoms. We remark that localized 
Fg7 orbitals with symmetry have linear coupling with JT 
phonons with Eg symmetry, since the symmetric representa- 
tion of eu X eu includes Eg. Then, H^-ph is given by 

-ffeph = 5(Q2r, + + (P| + P|)/2 

(15) 

+ {uj^/2){Ql + Ql) + HQl - IQlQ^), 

where g is the electron-phonon coupling constant, Qi and Qs 
are normal coordinates for [p? — y^)- and {iz^ — r^)-type JT 
phonons, respectively, P2 and P3 are corresponding canoni- 
cal momenta, T^= Y^aiflafha + fLfe^a), T^= J^aiflaha " 
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fl^fha), w is the frequency of local JT phonons, and b indi- 
cates the cubic anharmonicity. Note that the reduced mass for 
JT modes is set as unity. Here we introduce non-dimensional 
electron-phonon coupling constant a and the anharmonic en- 
ergy (3 as a=g'^/{2uj^) and f3=b/ {2u))^/'^, respectively. 

3.2 Multipole Susceptibility 

In order to clarify the magnetic properties at low temper- 
atures, we usually discuss the magnetic susceptibility, but 
in more general, it is necessary to consider the susceptibil- 
ity of multipole moments such as dipole, quadrupole, and 
octupole.^^ The multipole operator is given in the second- 
quantized form as 



(16) 



where X denotes the symbol of multipole with the synraietry 
of and 7 indicates a set of indices for the irreducible rep- 
resentation. For j=5/2, we can define multipole operators up 
to rank 5 in general, but we are interested in multipole prop- 
erties from the Fg quartet. Thus, we concentrate on multipole 
moments up to rank 3. 

Now let us show explicit forms for multipole operators fol- 
lowing Ref. 63. See also Ref. 62. As for dipole moments with 
symmetry, we express the operators as 

J^ux ~ Jxi Jiuy ~ Jyi Jiuz ~ ^zi (17) 

where J^, Jy, and are three angular momentum oper- 
ators for j=5/2, respectively. Concerning quadrupole mo- 
ments, they are classified into T^g and T^g. We express Fsg 
quadrupole operators as 

Osgu = (2J^ - - J^)/2, 

03g. = V3(Jj-J,2)/2. ^''^ 

For the quadrupole, we have the three operators 

Osg^ = \/3 JyJz/^., 

O^gr, = V3XX/2, (19) 

OsgC = JxJy/2, 

where the bar denotes the operation of taking all possible per- 
mutations in terms of cartesian components. 

Regarding octupole moments, there are three types as r2u, 
r4u, and Fsu. We express the octupole as 

T2u = \/T5XV^/6. (20) 
For the octupole, we express the operators as 



TauX — {"^Jx JxJy JxJz)/'^1 

while the octupole operators are given by 



(21) 



Tsuz 



15(J^ - J. J|)/6, 
15(J, J| - J, Jj)/6, 



(22) 



^{JzJl - JzJ^)/Q. 

We redefine the multipole moments so as to satisfy the or- 
thonormal condition Tr(X^X-y/)=^^-y/,^"' where 5^^i is the 
Kronecker's delta. 

In principle, the multipole susceptibility can be evaluated 
in the linear response theory, but we should note that the 



multipole moments belonging to the same symmetry can be 
mixed in general. In order to determine the coefficient of such 
mixed multipole moment, it is necessary to find the optimized 
multipole state which maximizes the susceptibility. Namely, 
we define the multipole operator as 



(23) 



where is determined by the eigenstate with the maximum 
eigenvalue of the susceptibility matrix, given by 



1 \ ^ e 

X77' = ;^2^- 



--E„/T 



-Err, IT 



E„ 



En 



{n\X^\m){m\Xy\n), (24) 



where is the eigenenergy for the n-th eigenstate |n) and Z 
is the partition function given by ■2^=X^„ e~^"/^. 

3.3 Method 

In order to evaluate multipole susceptibihties, here we em- 
ploy a numerical renormaUzation group (NRG) method,*^' 
in which momentum space is logarithmically discretized to 
include efficiently the conduction electrons near the Fermi en- 
ergy and the conduction electron states are characterized by 
"shell" labeled by N. The shell of N=0 denotes an impurity 
site described by the local Hamiltonian. The Hamiltonian is 
transformed into the recursion form as 

Hn+1 = VAHn +tN^i 

^iVcr^^+i'^ + Cjy_|_j^cjvct), (25) 

where A is a parameter for logarithmic discretization, cjva 
denotes the annihilation operator of conduction electron in the 
A'^-shell, and t indicates "hopping" of electron between N- 
and {N + l)-shells, expressed by 

(l + A-i)(l- A-^-i) 



In 



A 



-2N-3 



(26) 



2^/(1 -A-2^-i)(l 
The initial term Hq is given by 

Ho = A-l/2[iJeff + ifeph + V V{clj,^ + flco.)]. (27) 



The component of multipole susceptibility eq. (24) is evalu- 
ated by using the renormalized state. The free energy F for / 
electron is evaluated by 



F = -T lim 



InTre--^"/"^ - InTre"-^"/^ 



(28) 



Note that the temperature T is defined as T=A-(^-i)/2 
in the NRG calculation. The entropy Simp is obtained by 
Simp=—dF/dT and the specific heat Cimp is evaluated by 
Cinip=—Td'^F/dT^. In this paper, A is set as 5 and we keep 
3000 low-energy states for each renormaUzation step. The 
phonon basis for each JT mode is truncated at a finite number 
A^ph, which is set as A^ph=20. 

4. Results 

4. 1 Without Jahn-Teller phonons 

First let us consider the case without .ffeph- In Fig- 2(a), 
we show the multipole susceptibility x for U=\ and A=0.1 
with the ground state. At low temperatures, we find the 4u 
magnetic state and other multipoles are suppressed. Note that 
the 4u magnetic moment is expressed as 

Miua = PaJiua + (laTiua, (29) 



6 J. Phys. Soc. Jpn. 



Full Paper 



Takashi HOTTA 



(a) n=5 W=-6x10 V=0.05 
>.=0.1 x=0.4 

u=i .0 y=o.3 




(c) n=5 


W=-6x1 0"* 


V=0.05 


>.=0.1 


x=0.4 




U=2.5 


y=0.3 




2u 












■ 5u 














(d) 






log6- 


Jqg_4__ 


entropy 
V 








specific heat ^ ♦ ^ 



lo'^d' 10° Id' 10" id' Id' 10 10 10 id°id id id id" id' id" id' id' id' 

T T 



Fig. 2. (a) Tx and (b) Simp and C'imp vs. temperature for A=0.1 and 
(7=1.0. (c) Tx and (d) Simp and Cimp vs. temperature for A=0.1 and 
(7=2.5. 



for a=x, y, and z, since the moments belonging in the same 
symmetry are mixed in general.*''* For the present parameters, 
we find pa=Q326 and qa=— 0-946. As shown in Fig. 2(b), 
we find a residual entropy log 2 of doublet, but it is not 
the true ground state. For instance, when we increase the hy- 
bridization V, the 4u magnetic moment is screened even in 
the present temperature range. In any case, for the case of 
ground state, there are no significant higher-order multipoles. 
We find the standard Kondo physics, when the extra coupling 
to JT phonons is not considered. 

On the other hand, the low-temperature state is drastically 
changed, when the ground state is Fg^ quartet, as shown in 
Fig. 2(c). After the diagonalization of the susceptibility ma- 
trix, we obtain that the dominant multipole is purely given by 
2u octupole. The second and third multipole states have 3g 
quadrupole and 4u magnetic moments, respectively. The 4u 
moment is expressed by eq. (29) with pa=0. 560 and (7^=0.828. 
The fourth multipole state is given by 5u octupole, but in the 
case of n=5, 4u moment is not mixed with 5u, in contrast to 
the case of n=3P As explained in the previous section for 
the CEF states of n=2, 4u and 5u states in Oh symmetry are 
mixed in Th- Thus, the mixture of 4u and 5u moments is one 
of characteristic points of filled skutterudites with Th sym- 
metry. However, this mixing occurs due to the CEF term of 
Bq, which does not appear in the J=5/2 space for n=l and 5. 
Thus, in Sm-based filled skutterudites, the mixing of 4u and 
5u moments does not occur within an electronic model. In 
Fig. 2(c), we also find another 4u state and 5g quadrupole. 

In the case of the Fgy ground state, since this quartet carries 
all the multipole moments up to rank 3, we find the signifi- 
cant values for all moments. As shown in Fig. 2(d), we find a 
residual entropy of log 4, corresponding to the localized Fg^ 
quartet, since we have considered the hybridization between 
Ou conduction band and F^ state. In actuality, there should 
exist a finite hybridization between Cu conduction bands and 
F^7 states, even if the value is not large compared with that 
between conduction and F^ electrons. Thus, the entropy 
of log 4 should be eventually released at low temperatures in 



the actual situation. 

4.2 With Jahn-Teller phonons 

Next we include the effect of dynamical JT phonons, but 
before proceeding to the numerical results, let us briefly re- 
view the effect of JT phonons in the adiabatic approximation 
by following Ref. 28. Note that in actuality, the potential is not 
static, but it dynamically changes to follow the electron mo- 
tion. For /3=0, the potential is continuously degenerate along 
the circle of the bottom of the Mexican-hat potential. Thus, 
we obtain double degeneracy in the vibronic state concerning 
the rotational JT modes along clockwise and anti-clockwise 
directions. When a temperature becomes lower than a char- 
acteristic energy T*, which is related to a time scale to turn 
the direction of rotational JT modes, the entropy log 2 should 
be eventually released, leading to Kondo-like behavior, since 
the specific rotational direction disappears due to the average 
over the long enough time. 

In Figs. 3(a) and 3(b), we show multipole susceptibilities, 
entropy, and specific heat for cj=0.3, q;=0.5, and /3=0.0. Note 
that here we suppress the cubic anharmonicity. We find that all 
the multipole moments are still active, although the dominant 
multipole is changed from 2u to 4u. In this sense, the results 
do not seem to be qualitatively changed between Figs. 2(c) 
and 3(a), in spite of the JT active situation. We find a resid- 
ual entropy of log 4 in Fig. 3(b). In Fig. 3(c), we show the 
temperature dependence of average displacements. For the 
whole temperature range, we find finite values of v(Qf) and 
-\/ (Q§), while {Q2)={Qj)=0- Namely, JT vibrations occur 
around the origin without displacements. As the temperature 
is decreased, such vibrations are expected to be suppressed, 
but in the temperature range, we cannot find it. Rather, due to 
the isotropic JT vibrations, orbital fluctuations are still active 
in addition to magnetic fluctuations, leading to active multi- 
pole fluctuations. 

As explained in Ref. 28, when we include the effect of an- 
harmonicity, three potential minima appear in the bottom of 
the JT potential in the adiabatic approximation. Since the ro- 
tational mode should be changed to the quantum tunneling 
among three potential minima at low temperatures, the fre- 
quency is effectively reduced in the factor of e~^^l^ , where 
6E is the potential barrier. Then, we expect to observe the 
quasi-Kondo behavior even in the present temperature range, 
when we include the effect of cubic anharmonicity. 

In Figs. 3(d) and 3(e), we show multipole susceptibilities, 
entropy, and specific heat for aj=0.3, q:=0.5, and /3=— 0.006. 
For r>10~^, multipole susceptibilities are quite similar to 
Fig. 3(a). However, when temperature is decreased, the mixed 
multipole state with 4u magnetic and 5u octupole moments 
becomes dominant. This mixed moment is expressed as 



(30) 



for a=x and y. In the low-temperature region, we obtain 
Pa=0.761, qa=0.428, and ra=-0.488. Due to the effect of JT 
phonons, we find a significant contribution of 5u octupole. 
It is emphasized that the 5u component first appears in the 
present case for n=5. The secondary moment is given by 3gii 
quadrupole and the third state is given by 



= PzJiu 



(31) 
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tipole susceptibility, specific heat, entropy, average displace- 
ment for the ground state, we find the results quite similar 
to Figs. 3. In the present parameters, the ground state charac- 
ter is masked due to the effect of JT phonons. If we increase 
the CEF excitation energy, we could find different results from 
those of Fig. 3 for the ground state, but it is not so impor- 
tant in the context of Sm-based filled skutterudites with the 
CEF excitation in the order of lOK. 



(b) 
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Fig. 3. (a) Tx, (b) Si^p and d^p, and (c) (Q,) and ^JiQf) (i=2 and 3) 
vs. temperature for A=0.1, (7=2.5, and /3=0. (d) Tx, (e) Simp and Cimp, 
and (f) {Qi) and ^J{Qf) (i=2 and 3) vs. temperature for A=0.2, (7=1.0, 
and /3=-0.006. 



where pz=0. 67, q^=— 0.739, and rz=Q-0. Note that is found 
to be zero within the numerical precision, but it is difficult to 
conclude that it vanishes exactly. 

In Fig. 3(e), around at r=10^^, we find a peak in the spe- 
cific heat, since an entropy of log 2 is released. As mentioned 
above, this is considered to be quasi-Kondo behavior, origi- 
nating from the suppression of the rotational mode of dynam- 
ical JT phonons. In this case, an entropy of log 2 concern- 
ing orbital degree of freedom coupled with JT phonons is re- 
leased, while there still remains spin degree of freedom in the 
localized T^^ quartet. In fact, at low temperatures, magnetic 
susceptibility becomes dominant. 

In Fig. 3(f), we show the temperature dependence of av- 
erage displacements. For r>10^^, we find y/{Qf)^0 and 
V (QDt^O' suggesting that both Q2 and Q3 modes are active. 
This is consistent with the finite values of susceptibilities for 
Osgu and O^gv ■ Note that the displacement is not considered 
to occur, since {Q2)={Q3)=0- In the low-temperature region, 
after the release of an entropy log 2 of the rotational degree of 
freedom, we find y(g|y=(Q2)=0, while ^/J^=\(Q-i)\^0, 
indicating that only Qs-type JT vibration is active with fi- 
nite displacement. This is also consistent with the result that 
X3gu remains at low temperatures, since the vibration mode is 
fixed as Qa-type, after the quasi-Kondo phenomenon occurs. 
Namely, the JT vibrations become anisotropic for T<T*, 
while they are isotropic for T>T*, as mentioned in Ref. 28. 

For the case with JT phonons, we have shown only the re- 
sults for the Fgy ground state. When we evaluate the mul- 



5. Discussion and Summary 

In this paper, we have examined the multipole state of Sm- 
based filled skutterudites on the basis of the multiorbital An- 
derson model with the use of the NRG method. As mentioned 
in the introduction, a possibility of octupole ordering has been 
discussed experimentally in SmRu4Pi2.^ " In general, there 
are two possibilities as 2u and 5u octupoles. The 2u octupole 
has been found to be dominant in the electronic model, when 
we do not include the coupling between / electrons and JT 
phonons. On the other hand, the 5u octupole becomes signif- 
icant in the mixed state with 4u magnetic moment, when we 
consider the coupling to JT phonons. 

In the experiment, Yoshizawa et al. have suggested the 
5u octupole from the elastic constant measurement.'" The 
isotropic 2u octupole seems to contradict with experimental 
results. Note here that the pure electronic model cannot stabi- 
lize the 5u octupole, since there is no mixing between 4u and 
5u moments for the case of n=5 with the ground state multi- 
plet characterized by J =5 12. To obtain the mixed multipoles 
of 4u magnetic and 5u octupole, it is essentially important 
to consider the electron-phonon coupling. It is not yet clari- 
fied that the dominant phonon mode in filled skutterudites is 
Jahn-Teller type with Eg symmetry, but in any case, the effect 
of rattling seems to play a crucial role for the appearance of 
the octupole state. It may be interesting to carry out an ex- 
periment to detect a coupling between rattling and /-electron 
state in SmRu4Pi2. 

Concerning the mechanism of magnetically robust heavy- 
fermion phenomena observed in SmOs4Sbi2,'^ a potential 
role of phonons has been pointed out from the viewpoint of 
the Kondo effect with non-magnetic origin.^'' In this context, 
the quasi-Kondo behavior due to the dynamical JT phonons 
may be a possible candidate to understand magnetically ro- 
bust heavy-fermion phenomenon. Further investigations are 
required, in particular, concerning an experimental method to 
detect the Kondo-like behavior due to Jahn-Teller phonons. It 
is one of important future issues. 

In summary, we have discussed the multipole state for the 
case of n=5 by analyzing the multipole Anderson model with 
the use of the NRG technique. When we do not consider the 
coupling between JT phonons and / electrons in F^y quartet, 
we have found that the dominant multipole moment is 2u oc- 
tupole. When the coupling with JT phonons is switched and 
the cubic anharmonicity is included, the dominant multipole 
state is characterized by the mixture of 4u magnetic and 5u oc- 
tupole moments, with significant fluctuations of quadrupole. 
We have also observed the quasi-Kondo behavior due to the 
entropy release concerning the rotational JT mode. 
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